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Abstract. We study parton suppression effects in heavy-ion collisions within the Parton Quenching Model 
(PQM). After a brief summary of the main features of the model, we present comparisons of calculations 
for the nuclear modification and the away-side suppression factor to data in Au+Au and Cu+Cu collisions 
at ^/s NN = 200 GeV. We discuss properties of light hadron probes and their sensitivity to the medium 
density within the PQM Monte Carlo framework. 



1 Introduction 

One of the very early and very exciting findings from the 
experiments running at the Relativistic Heavy Ion Col- 
lider (RHIC) was the observation of apparent jet suppres- 
sion. At the top RHIC energy, ^/s NN = 200 GeV, the mid- 
rapidity yield of high transverse momentum leading par- 
ticles in Au+Au collisions is about a factor of five lower 
than expected from the measurements in p+p collisions at 
the same energy [J, (2|. Similarly, jet-like correlations on 
the azimuthally-opposite ('away') side of a high-p t trigger 
particle are suppressed by a factor of four to five, while 
the near-side correlation strength is almost unchanged fj| . 
The observed suppression persists also at lower center-of- 
mass energies H, as well as in Cu+Cu collisions at 62.4 
and 200 GeV [5|. The absence of these effects in d+Au 
collisions at ^/s NN = 200 GeV Q strongly supports the 
picture of partonic energy loss, where energetic partons, 
produced in initial hard scattering processes, lose energy 
as a consequence of the final-state interaction with the 
dense partonic matter created in nucleus-nucleus colli- 
sions. The dominant contribution to the energy loss is 
believed to originate from medium-induced gluon radia- 
tion (see Ref. [7| and references therein). Recent calcula- 
tions point out, however, that the collisional contribution 
to the energy loss might not be negligible [§J. Still, strong 
interest in these probes arises mainly from the fact that 
modification of their properties due to interaction with 
the medium provides access to fundamental properties of 
the created matter, such as its density and nature [§[]. 

A quite simple model that includes final-state gluon 
radiation is the Parton Quenching Model (PQM) [l0|. 
It combines the pQCD BDMPS-Z-SW framework for the 
probabilistic calculation of parton energy loss in extended 
partonic matter of given size and density jTlj ] with a re- 
alistic description of the collision overlap geometry in a 
static medium. High-energy partons (and parton pairs) 



are treated on an event-by-event basis using Monte Carlo 
techniques. The model has one single free parameter that 
sets the scale of the medium transport coefficient q, the 
average transverse momentum squared transferred to the 
hard parton per unit path length, and, thus, the scale of 
the energy loss. 

In these proceedings, we will, after a short introduction 
to the PQM model, compare the suppression phenomena 
introduced above to the calculations obtained with PQM. 
In particular, we will discuss to what extent light hadronic 
probes are sensitive to the dense matter formed at RHIC. 



2 The PQM model 

The Monte Carlo calculation of the unquenched and quen- 
ched transverse momentum spectra in PQM consists of 
four main steps: 1) Determination of a parton type and 
its transverse momentum according to PYTHIA (LO) par- 
ton distribution functions. 2) Determination of its parton- 
production point in the transverse plane according to the 
nuclear density profile (Glauber) and evaluation of path 
length and transport coefficient seen by the produced par- 
ton using two line integrals weighted with p co i\. 3) Cal- 
culation of the energy loss using constrained quenching 
weights to extrapolate from the eikonal approach used 
in the BDMPS-SW framework to finite parton energies. 
Two types of constraints are constructed to estimate the 
systematic uncertainty of the approach: In the reweightcd 
case the energy-loss distribution is simply normalized with- 
in the kinematically allowed regime, whereas in the non- 
reweighted case the fraction of the distribution larger than 
the energy of the parton contributes to the probability of 
maximum energy loss. 4) Finally, independent fragmenta- 
tion is applied to the quenched and original parton. Back- 
to-back parton pairs initially consist of a pair of partons 
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with the same pt at the same production point, but with 
opposite-side emission angle. 

Currently, the model is restricted to mid-rapidity, and, 
for simplicity, we ignore initial-state effects. The single 
parameter of the model (k) is fixed to set the scale of the 
energy loss by fitting to data from 0-10% central Au+Au 
collisions at ^/s NN = 200 GeV. Once the scale is fixed we 
implicitly vary the medium density by its dependence on 
the centrality as given by Glauber. 

Details on the quenching procedure and its application 
to high-pt data can be found in Ref . [Kj . 



3 Suppression of leading particles 

Suppression of leading particles is usually quantified via 
the nuclear modification factor, 

(A'coii) d 2 iVpp/dptd?7 

the ratio of the yield of light hadrons in nucleus-nucleus 
over proton-proton yield scaled by the number of binary 
collisions in a given centrality class. The ratio is normal- 
ized so that, if no final-state effects were present, it would 
be close to one. Indeed, this scaling is observed in the 
measurement of dire ct p hoton yield in Au+Au collisions 
at ^/i^N = 200 GeV pjj. 

However, in Au+Au collisions at the same energy, Raa 
for light hadrons at mid-rapidity is found to decrease from 
peripheral (Raa — 1) to central events (Raa — 0.2), for 
Pt ~ 5 GeV (see Fig. UJ). In this high-p t region, Raa is in- 
dependent of the particle type and rather independent of 




Fig. 1. 7?aa(p0 in Au+Au at = 200 GeV for differ- 

ent centralities. Data are PHENIX charged hadrons and tt° [3] 
and STAR charged hadrons [1] with combined statistical and 
Pt-dependent systematic errors (bars on the data points) and 
Pt-independent systematic errors (bars at -Raa = !.)■ The gray 
band is the original PQM calculation from Ref. [lOj using re- 
weighted and non-reweighted quenching weights. 



pt. In Fig. [TJ the data is compared to the original calcu- 
lation of PQM, where the single parameter of the model 
was adjusted to match the suppression measured in 0-10% 
central Au+Au collisions at ^/s NN = 200 GeV. For central 
collisions, this corresponds to an average transport coeffi- 
cient of about (q) = 14 GeV 2 /fm, which decreases with de- 
creasing centrality to essentially zero for the most periph- 
eral collisions. The average is taken over all produced hard 
partons and given in the equivalent static scenario [jT| . 

In Fig. \2\ we show i?AA(j>t) for neutral pions in 0- 
10% central Au+Au collisions at ^/s NN = 200 GeV for pre- 
liminary PHENIX data [T3j, which extends beyond the 
previous data by almost a factor of two in p t . The new 
data is compared to the original PQM calculation for 
(q) = 14 GeV 2 /fm (simply extended to larger p t ) in the re- 
weighted and non-reweighted approximation of the quench- 
ing weights for finite (small) parton energies. For a given 
medium density, the reweighted values generally produce 
a weaker energy loss than the non-reweighted ones. How- 
ever, in the reweighted case, only a rather small fraction 
of the full quenching weights is used. This results in the 
unphysical behaviour that, beyond a certain medium den- 
sity, the probability of 'zero energy loss' increases for in- 
creasing densities. Therefore, we show in addition, results 
for lower densities, for (q) = 4 and 7 GeV 2 /fm, in the 
non-reweighted case. We furthermore compare the realis- 
tic calculations with calculations where we either fix the 
geometry (L — 6 fm and q = 1 GeV 2 /fm) or the relative 
energy loss (AE/E = 0.25). While the fixed geometry case 
leads to a more strongly rising Raa with increasing pt, the 
fixed relative energy loss approximation over a very wide 
range of p t describes the data. 

Similar conclusions can be made for the comparison 
with data from Cu+Cu interactions at RHIC. In Fig. [3] we 
show i?AA(pt) for neutral pions in 0-10% central Cu+Cu 
collisions at ^/s NN = 200 GeV for preliminary PHENIX 
data [l3| . The data are compared to the PQM prediction 
(made before QM'05 [l4j) for (q) = 9 GcV 2 /fm for the 
reweighted and non-reweighted case. The value of (q) = 
9 GeV 2 /fin relies on the proportionality of the transport 
coefficient to the initial volume-density of gluons [HJ and 
on the predictions of the saturation model [16| , as outlined 
in Ref. [Io| • The prediction seems to slightly overestimate 
the suppression in the smaller Cu+Cu system. We there- 
fore also present results for lower densities, for (q) = 3 and 
5.5 GeV 2 /fm, in the non-reweighted case. As before, the 
calculations including the nuclear geometry are compared 
to calculations with either fixed geometry (L — 4 fm and 
q = 1 GeV 2 /fm) or relative energy loss (AE/E = 0.15). 

In Fig. 3} we show calculations of the centrality depen- 
dence of the nuclear modification factor for p t > 7 GeV in 
Cu+Cu and Au+Au collisions at ^/s NN = 200 GeV com- 
pared to the preliminary PHENIX data [l3j. The data 
seem to favour values for the transport coefficient of (q) = 
7-14 for Au+Au and 3 - 5.5 GeV 2 /fm for Cu+Cu colli- 
sions. 



Constantin Loizides: High transverse momentum suppression and surface effects in Cu+Cu and Au+Au collisions 3 



< 
< 
cr 



0-10% Au+Au, 200 GeV 



1-2 — 

t 

0.8 - 



0.6 - 



PHENIX (prel., QM05)jl° 

<q> = 4, 7 and 14 GeV 2 /fm, nrw 

<q>=14GeV 2 /fm, rw 

fixed L = 6 fm, q = 1 GeV 2 /fm, nrw 

fixed AE/E =0.25 




20 25 

P, [GeV] 

Fig. 2. i?AA(pt) for neutral pions in 0-10% central Au+Au 
collisions at ^/s NN = 200 GeV (preliminary PHENIX data) 
with combined statistical and pt-dependent systematic er- 
rors, as well as p t -independent systematic errors (bar at 
Raa = 1) 0]. The PQM curves for (q) = 14 GeV 2 /fm are the 
original PQM results (extended to larger pt) for the reweighted 
and non-reweighted case from Ref. [101 ] . In addition, results for 
(q) = 4 and 7 GeV 2 /fm for the non-reweighted case are shown, 
as well as calculations for fixed L = 6 fm and q = 1 GeV 2 /fm, 
or for fixed relative energy loss of AE/E — 0.25. 
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Fig. 3. Raa(jPi) for neutral pions in 0-10% central Cu+Cu col- 
lisions at ^/s NN = 200 GeV (preliminary PHENIX data) with 
combined statistical and pt-dependent systematic errors, as 
well as pt-independent systematic errors (bar at Raa = 1) fl3| . 
The PQM curves for (q) = 9 GeV 2 /fm are the original PQM 
predictions for the reweighted and non-reweighted case from 
Ref. flQj . In addition, results for (q) — 3 and 5.5 GeV 2 /fm for 
the non-reweighted case are shown, as well as calculations for 
fixed L — Aim and q = 1 GeV 2 /fm, or for fixed relative energy 
loss of AE/E = 0.15. 
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Fig. 4. Centrality dependence of 7?AA(pt > 7 GeV) for neutral 
pions in Cu+Cu and Au+Au collisions at -^/s NN = 200 GeV 
(preliminary PHENIX data) with combined statistical and sys- 
tematic errors fl3l ]. The PQM calculation is shown for (q) = 4, 
7 and 14 GeV 2 /fm (Au+Au), as well as for (q) = 3, 5.5 and 
9 GeV 2 /fm (Cu+Cu) in the non-reweighted case. 
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Fig. 5. Centrality dependence of /^A ay m Au+Au at ^/s NN = 
200 GeV for two data sets. The data for the low p t cuts of 
4 < p' rig < 6 GeV and 2 GeV < p t assoc < pT ig with statistical 
(bars) and systematic (ticks) errors are from Ref. (3) . The data 
for the high p t cuts of 8 < pf ig < 15 GeV and 6 GeV < pt ssoc < 
p* rlg with statistical (bars) errors are from Ref. [2C| (wrt d+Au 
instead of pp). PQM results in the non-reweighted case for 
Au+Au are for (q) = 4, 7 and 14 GeV 2 /fm, and for Cu+Cu 
for (q) — 3, 5.5 and 9 GeV 2 /fm (only high trigger cuts). 



4 Suppression of jet-like correlations 

Within the PQM framework, we consider the suppression 
of back-to-back jet-like correlations by simulation of back- 
to-back pairs of partons in a simple LO parton model (no 
intrinsic k t ). The magnitude of the suppression is usually 
quantified by the factor Iaa^ = ^aa Y /^pp ay > where the 
di-hadron correlation strength, -Dpp^.A)' ^ or an assocrai, ed 
hadron, /12, with p+,,2 in the opposite azimuthal direction 



from a trigger hadron, hi, with p t ,i, is integrated over the 
considered trigger- and associated- p t intervals [13] , 

^MAA)=/!f X d PM /27dp t , 2 (2) 

r d3 %p(l A )/dpt,idp t , 2 d^ 

./away side df7 pp(AA) /^Pt,! 
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As was found for the nuclear modification factor in 
Au+Au collisions at ^/s NN = 200 GeV, Ij^ y for 4 < 

pl lig < 6 GeV and 2 GeV < p assoc < p£ rig is found to 
decrease with increasing centrality, down to about 0.2-0.3 
for the most central events (see Fig. [5). In the figure, we 



also show data taken from Ref. [2(jJ for higher p t 

cuts of 8 < p" s < 15 GeV, where we normalize the away- 
side yields measured in Au+Au to the yields measured 
in d+Au collisions. Both data sets are compared to PQM 
calculations (using p+p as the reference in both cases) for 
the non-reweighted case with (q) = 4, 7 and 14 GeV 2 /fm 
for Au+Au collisions. The data for the higher trigger cuts 
favour smaller medium densities, while the uncertainties 
on the data for the lower cuts make it difficult to draw 
strong conclusions. For completeness, we also present pre- 
dictions for the Cu+Cu system with (q) = 3, 5.5 and 
9 GeV 2 /fm (only high trigger cuts). 



5 Sensitivity of light hadronic probes 
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Fig. 6. Raa at 10 GeV and I^ y for trigger cuts of 8 < pf ig < 
15 GeV and 6 GeV < pf soc < pf ig as a function of (q) in 0- 
5% central Au+Au and Cu+Cu collisions at ^/s NN = 200 GeV. 
The calculations are in the non-reweighted approximation. 



Another complication arises from the fact that probes 
based on leading-particle analyzes arc affected by sev- 
eral biases. The biases result primarily from the steeply- 
falling underlying production cross sections and the emis- 
sion from regions close to the surface. The latter effect 
dominates for large medium densities. This is illustrated 
in Fig. \6\ which shows the behaviour of Raa (at 10 GeV) 
and Iaa Y (f° r trigger cuts of 8 < p\ lls < 15 GeV and 
6 GeV < p assoc < pi" 6 ) as a function of the medium den- 
sity, expressed as (q), in Au+Au and Cu+Cu collisions 
at -y/s NX = 200 GeV. Clearly, beyond a certain medium 
density, the numerical values for the ratios saturate at a 
non-zero value. The fundamental reason is that the proba- 
bility of no medium- induced gluon radiation, P(AE = 0), 
for a medium of finite size and finite density is not zero. 
The 'no radiation' contribution to the spectra is dramat- 
ically enhanced when realistic nuclear path-length distri- 
butions and density profiles are taken into account, since 
P(AE = 0) is decreasing as a function of qL 3 , giving sig- 
nificant weight to partons that 'feel' lower values of qL 3 . 
We have reported in Ref. [k| that the Raa and /AA ay data 
(for low trigger cuts) can be described by taking into ac- 
count only P(AE = 0) and 1 - P{AE = 0) (see also [l8j). 

Recently, direct measurements of dijets in heavy-ion 
collisions have been performed by the STAR collabora- 
tion [20| . It is found that, while the relative yields of 
the hadron-triggcrcd fragmentation function relative to 
d+Au, the integrand in eq. (2), are suppressed, the shape 
is not modified, even in the most central Au+Au collisions 
at 200 GeV. Such a scenario naturally follows from the 
trigger bias since, due to the cuts on near-side and away- 
side particle p t , those dijet pairs that escaped the collision 
region without losing a significant fraction of their initial 
energy are preferentially selected. 

For the PQM framework, the situation is illustrated 
in Fig. [7J which shows production points and emission di- 
rection for surviving back-to-back pairs in the transverse 
plane for 0-10% central Au+Au collisions. Each parton 
pair shown yields a hadron pair within 8 < p t rlg < 15 GeV 



and 6 GeV < p assoc < pf ig . The chosen set of trans- 
port coefficients are the same as previously used to de- 
scribe the -?AA ay data, (q) = 4, 7 and 14 GeV 2 /fm, in 
the non-reweighted case. In panels of Fig. \7\ the center 
of any line indicates the production point of a parton 
pair, and the two partons emerge along the line in oppo- 
site direction. The color of each line either indicates the 
medium density (left) or the relative energy loss (right 
panels) of the away-side parton, which is defined as the 
parton fragmenting into the associated hadron. The left 
panels of Fig. \7] clearly illustrates that the largest medium 
densities are encountered by away-side partons that pass 
through the central zone of the collision, which is ex- 
pected for a static scenario. The core is surrounded by 
dijets that escape tangentially with respect to the cen- 
tral region Ref. 19j- In the case of an expanding medium, 
the interpretation of the underlying trigger bias remains 
the same: surviving dijets are always selected from regions 
(and times) of the collision evolution where they suffered 
the least amount of interactions with the medium. The 
right panels of Fig. \7\ which color-code the relative en- 
ergy loss of the away-side parton, illustrate this point from 
a slightly different perspective. It is apparent that in all 
three cases, the energy loss is close to zero for a large 
fraction of the dijets, even for the core region. This is 
due to the probability of zero-energy loss for the away- 
side parton, p(AE away = 0), which for the core is already 
about 0.05, even for the (q) = 14 GeV 2 /fm case. As shown 
in Fig. |8(a)[ P(AE away = 0) quickly reaches values be- 
tween 0.1 and 0.3 for the d^4 fm, where d = \J x\ + iJq 
is the distance of the production point to the collision 
center. A large fraction of these parton pairs are initially 
produced, see Fig. 8(b) (dashed line, no medium). The 
same figure illustrates the shift of the production point 
distribution to larger distances with increasing medium 
densities, when energy loss is included. The mean values 
change from 3.2 fm in vacuum to 3.5, 3.7 and 4.1 fm for 
(q) = 4, 7 and 14 GeV 2 /fm, respectively. In Fig. |8(c)[ 
we finally present the average AE/E versus d for near- 
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Fig. 7. (Color online) Production points and emission direction for surviving back-to-back parton pairs (yielding hadron 
pairs within 8 < pf lg < 15 GeV and 6 GeV < p^ ssoc < p* rlg ) in the transverse plane for (q) = 4 (bottom), 7 (middle) and 
14 GeV 2 /fm (top panels) for 0-10% central Au+Au collisions. The line color represents the medium density (relative energy 
loss) of the away-side parton in the left (right) panel. 
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side and away-side partons. For distances of d > 4 fm, that 
contribute significantly to the near- and away-side yields, 
0.1 < AE/E < 0.25 for the three densities, and the average 
relative energy loss is nearly the same on both sides. 



6 Summary 

Jet quenching effects at the top RHIC energy are discussed 
within the Parton Quenching Model that includes a prob- 
abilistic treatment of the BDMPS quenching weights and 
a Glauber-based implementation of the collision geometry 
The available high-p t data for Raa, Iaa Y an( ^ their cen- 
trality dependence constrain the extracted medium den- 
sity to about 3 ;$ (q) ^9 GeV 2 /fm for central Cu+Cu and 
4< (q) < 14 GeV 2 /fm for central Au+Au collisions. Our 
analysis suggests that particle production in central colli- 
sions is 'surface' dominated, not only for single hadrons, 
but also for dijets. The properties of surviving dijets there- 
fore arc very similar to the vacuum properties. 
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Fig. 8. (a) The probability of zero-energy loss for the away- 
side parton versus the distance of the production point from 
the collision center, [(b)] the distribution of production distances 
from the collision center and |(c)| the average relative energy loss 
for near-side and away-side partons versus the distance of the 



